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Abstract 

A large- Nf expansion is used to compute the partially quenched chiral condensate of QCD in 
the microscopic finite- volume scaling region. 
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While the replica method has found widespread use in condensed matter physics as a means of gen- 
erating quenched averages, it has only rarely been applied in particle-physics contexts. Analytical 
approaches to quenching - so often used in lattice gauge theory simulations - have typically been 
based on the alternative supersymmetric method as applied to the effective chiral Lagrangian (see 
also ref. 0). In the past few years this supersymmetric method has been successful in deriving exact 
hnite-volume scaling laws for the partially quenched chiral condensate, and in computing analytically 
the smallest Dirac operator spectrum in gauge theories with spontaneous chiral symmetry breaking 
Q. Because the supersymmetric method relies on a subtle definition of the global integration super- 
manifold of the effective Lagrangian, it is nevertheless of interest if the same results can be arrived 
at in an entirely different manner. Steps in that direction were made very recently in ref. [|], where 
it was shown that small-mass and large-mass series expansions of the partially quenched chiral con- 
densate can be computed via the replica method. In the present paper we shall extend these results 
by showing that also a large-./Vy expansion is a suitable starting point for the replica method. The 
resulting series expansion for the partially quenched chiral condensate turns out to be a resummed 
version of the asymptotic large-mass expansion derived in ref. and thus more accurate in the 
appropriate regime. 

In the present context the replica method consists in adding to the SU(iV c > 3) QCD Lagrangian of 
Nf physical quark flavors N v additional "valence quarks", which, in the simplest case, are taken to 
be mass-degenerate. We denote the physical quark masses by m/, and the mass of the valence quarks 
by m v . We next consider the extended QCD partition function in a sector of topological charge v, 
which we for simplicity take to be non-negative (all physical results are in any case symmetric under 
v — > —v): 

(N f \ N f 

J] m) m^ v f [dA]„ det'(0 - m v ) N « det W " m f) e" 5 ™ [A] . (1) 
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The prime on the determinants indicate that they exclude zero modes, whose effect has already been 
taken into account by the prefactors. This is an (unnormalized) average of N v identical replicas of the 
fermionic partition function 



Z v = J dipdip exp J d 4 xip(i0 — m v )ip 



(2) 



Letting N v = just reproduces the original QCD partition function. But the theory extended with N v 
additional quark species can serve as a generating functional of partially quenched averages of ipip if 
one sets N v to zero after having performed the required differentiations. Assuming that J\f = Nf + N v 
is small enough for spontaneous chiral symmetry breaking to take place according to the conventional 
SU(A/ r ) x SU(AT) — > SU(AT) , the theory can be analyzed in the low-energy region by means of its effective 
chiral Lagrangian. For N v integer, this is completely straightforward, and no additional assumptions 
are required. In order to apply the replica method one must in addition rely on analyticity in the 
number of flavors M. This is the non-trivial step, and it is a priori far from obvious that it can be done 
at all. But as shown in ref. [|J|, the needed analytic continuation is possible in both small-mass and 
large-mass expansions. This is because the dimension of the coset group SU(A/") in these expansions 
appears only in simple group invariants that are at most rational functions of polynomials in AT. Except 
for the appearance of poles, such functions involving only polynomials of M can unambiguously be 
continued to non-integer A/". As in ref. Q, we shall restrict ourselves to the low-energy and finite- 
volume regime where V <C l/ m 7r i n ^ rie limit where V is sent to infinity. In this finite- volume scaling 
region there are also exact non-perturbative results from both the supersymmetric method || and 
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universal Random Matrix Theory formulas [|| with which to compare our results. The conventional 
infinite-volume chiral condensate is always normalized to S, independently of the physical number of 
flavors Nf. 



The (mass-dependent) partially quenched chiral condensate in the given finite four-volume V is defined 
by 

^iWl) s lim J_J_ lnZ CA0. (3) 

As indicated, in the above scaling region the masses only enter in the combinations \it = nifVE 
and /ij, = m v VT,. In lattice gauge theory simulations the implied finite- volume scaling has been 
checked with staggered fermions || and, very recently, also with lattice fermions sensitive to gauge 
field topology M. 



The effective finite-volume partition function is 



Z [M) 



U{N) 



dU (det U) u exp 



(4) 



where M is the J\f x TV matrix of rescaled masses: fj,f = m.f'EV for / = 1, . . . , Nf and /i„ = m v Y,V in 
the remaining entries. We shall here compute the partially quenched chiral condensate by means 
of a large- Aj expansion of this effective partition function ||. Because the effective action in eq. 
(||) involves A" terms, while the integration is over A" 2 variables, it is conventional to perform the 
1 /AA-expansion by considering an action rescaled by M: 



dU (det U) v exp \NTt{AU ] + UA*) 



(5) 



Letting A a denote the eigenvalues of A* A, the relation between these and the rescaled masses are thus 

At Nf = oo, the starting point of a large- Nf expansion, the effective partition function is known to 
describe two diffferent phases S. Defining 



I 



0fc 



-y 



(7) 



the transition between the two phases occurs at a\ = 2 for v = ||, and this remains unchanged for 
any finite v. We shall here focus on the phase with o\ < 2, which thus corresponds to large masses. 
We also note that 

^^^Ei. (8) 

a Pa 

and introduce in addition 

(9) 



Conventionally one defines the free energy by = {l/N)\nZ v , but we shall instead need 
J\fF^ = lnZ u (see e.g. eq. (|3|)). In the large- A/" expansion one most readily computes 



d\ a 



Ha dfJ, a 



(10) 



3 



which is particularly convenient for our purposes, since this is directly related to the partially quenched 
chiral condensate: 

E v (pv,{l*f}) Id 
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N v =0 



We thus need A^-fold degenerate eigenvalues (i v , and can trivially take the N v — ► limit (as was to 
be expected in an expansion around Nf = oo). 

In the sector of topological charge v = a large- M expansion of the effective partition function (||) 
has been worked out to high orders by Gross and Newman ||, and we can directly make use of their 
results. They find, in the phase with u\ < 2, 



16A/' 2 (2- ( j 1 )A: 
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(12) 



It is not a priory obvious that this large- M expansion is suitable for our purposes, since we are certainly 
only interested in the expansion for relatively small Nf. In fact, we shall even also seek the other 
extreme limit of Nf = 0, which just corresponds to the fully quenched case. However, in the phase 
with o\ < 2 the large-TV^ expansion ought to be connected with the saddle-point expansion, which in 
turn should be related to the large-mass expansion derived in ref. As we shall explicitly verify 
below, this is indeed the case. 



Using relations @ and (|T(]) we directly recover an expansion for the partially quenched chiral con- 
densate in a gauge field sector of vanishing topological charge: 



£o(/i„,{/V» 



N f 



i 



i 



+ 



189 



1 



+ 



9 



512(1 -ax) 



[1-ax) 128(1 -ffi) 4 
225 



(73 



r«"5 + ■ • 



1 

14 



u 6 



9 



1024(1 - ai) 6 
189 
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128(1 -cti) 3 512(1 -ax 
225 
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1024(1 - a x ) 



+ ... 



(13) 



Remarkably, all explicit factors of J\f have cancelled, and the above expansion is perfectly well-defined 
for any Nf, even all the way down to Nf = 0.Q The only requirement is that ax < 2, which translates 
into ax < 1 (in the Nf = case there is no such requirement, and one simply sets all void sums 
to zero in the above expansion). As expected, the resulting series is thus related to the large-mass 
expansion of ref. ||]. In fact, many curious regularities observed in the large-mass expansion of ref. 
Q now find their explanation: the above large-AA series is a non-trivial resummation of the ordinary 
large-mass expansion that highlights the special role played by the first inverse mass sum ax- In ref. 

1 Note that the limit N v = has already been taken here, so that Af — Nf at this point. 
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the large-mass series was for convenience truncated at an order (7th) taken to be the same both in 
the quenched mass \x v and the physical masses [if. Indeed, if one truncates the expansion of eq. ( |i~3| ) 
at the same order, one recovers exactly the large-mass series derived in ref. §. It is also in this form 
that one can most easily compare with the analytical results of both the supersymmetric method 
and Random Matrix Theory Q, and as already noted in ref. Q there is perfect agreement with both. 
However, only the large- Nj expansion has revealed that a resummation of the series is possible. 



The large-mass expansion of ref. @ made use of the fact that the v = effective partition function 
factors into two pieces, one originating directly from the leading-order saddle point solution, and a 
remaining function which is annihilated by a set of Virasoro generators. Unfortunately, it is not 
presently known how to extend this technique based on Virasoro constraints to the case of v ^ 0. But 
the large-TV expansion is not afflicted with this problem, and it can easily be carried through to the 
case v ^ 0. Also here much of the work has already been done for us. In particular, Brihaye and Rossi 
have derived the large- j\f expansion up to order 1//V 4 in the free energy. Using our notation their 
result for the free energy i?H = 1/vVln^ reads 



= 2£Ay 2 -^£M^ /2 + A 6 1/2 ] 
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This corresponds to 
p(y) - . 
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which in turn, introducing the rescaled sums (|9|) and the identification (|T 
quenched chiral condensate computed up to 0(l/Af 6 ): 



(15) 

leads to a partially 
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Of course, for u = we recover the leading orders of the expansion (13). We also note that the first 
terms in eq. (|l5|), those from the classical saddle point, are unaffacted by the presence of (det U) u in 
the integrand, as expected. If we expand the result ( |l6| ) in a\ we find 
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(17) 
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As a check, this agrees with the asymptotic expansion of the the fully quenched (Nf = 0) chiral 
condensate || |3[ 

" - = M« [I v (fiv)K v (ji v ) + I v+ i(ix v )K v -x(n v )] -\ , (18) 

as obtained from either the super symmetric method or from Random Matrix Theory. However, the 
standard asymptotic expansions of these expressions, using the individual asymptotic series of the 
modified Bessel function I n {x) and K n (x), miss the possibility of resummations as in the large- A/" 
expansions (|i~3| ) and (|l(^). It is interesting to note that the otherwise ubiquitous v/n v -texm. in the 
partially quenched chiral condensate, which normally [|l0| can be traced directly back to the factor 
from the zero modes in the original QCD partition function ([!]), is missing in the above expansion. It 
is due to a precise cancellation between this term and a (!///„) piece arising from the asymptotic 



expansions of the modified Bessel functions in (|18|), which thus, beyond the classical saddle-point 
terms, starts at C?(l//x^)^| 

To conclude, the large-AA expansion of the effective partition function of QCD in the finite-volume 
scaling regime V <C is a good starting point for an analytical computation of the partially 

quenched chiral condensate by means of the replica method. The resulting series is a resummation of 
the large-mass expansion derived by means of the same replica method in ref. M. In the u = sector 
we have explicitly checked this agreement, while for sectors with i/^Owe have confirmed in the Nf = 
case that the resulting series agrees with the asymptotic expansion of the partially quenched chiral 
condensate as it obtained by both the supersymmetric technique and Random Matrix Theory. The 
replica method is clearly a suitable alternative to these methods, at least in series expansions. To go 
beyond such expansions requires an analytical handle on "external-field" group integrals for unitary 
groups U(n) with n extended beyond integers, a problem that may appear forbiddingly difficult. 
Fortunately, less is actually needed: only the exact coefficient of the leading term in an e-expansion of 
the group \J(Nf + e), where Nf is integer. Finding this exact expression will allow one to establish the 
surprising relationship between such group integrals and those of the supergroup extensions discussed 
in ref. H. Other outstanding questions are how to use this method to derive spectral correlation 
functions of the Dirac operator Q , and how to understand why the correlations of the smallest Dirac 
eigenvalues can be described in terms of a theory extended with additional fermion species [jllj (a 
result that has a natural explanation in the supergroup formulation ||). Finally, the replica method 
may also find use as an alternative way to perform partially quenched chiral perturbation theory in 
general, i.e, beyond the finite-volume scaling regime considered here. 
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